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We present a search for the Higgs boson in the process qq —> ZH —> t^t~bb. The analysis uses 
an integrated luminosity of 1 fb^^ of pp collisions produced at y/s — 1.96 TeV and accumulated by 
the upgraded Collider Detector at Fermilab (CDF II). We employ artificial neural networks both to 
correct jets mismeasured in the calorimeter, and to distinguish the signal kinematic distributions 
from those of the background. We see no evidence for Higgs boson production, and set 95% C.L. 
upper limits on azH ■ B{H -^ bb), ranging from 1.5 pb to 1.2 pb for a Higgs boson mass [van) of 
110 to 150 GcV/c^ 
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The Higgs boson is the only particle predicted by the 
standard model (SM) of particle physics which has not 
yet been discovered. It is the physical manifestation 
of the mechanism which provides mass to fundamental 
particles [l|, Q. Direct searches have excluded the SM 
Higgs boson for masses mn < 114.4 GeV/c^ at the 95% 
C.L. [3|]. The Higgs boson mass is indirectly constrained 
from precise electroweak measurements to mn = 76_24 
GeV/c^ [J]. A number of extensions to the SM predict 
a SM-like Higgs boson, in particular Ref. [51] predicts a 
SM-like Higgs boson with 68% posterior probability to 
be between 115.4 and 120.4 GeV/c^. Only the Tevatron 
collider experiments are currently capable of extending 
the limits on a Higgs boson for niH > 114.4 GeV/c^. 

This Letter presents the first CDF II search for a Higgs 
boson in the process pp —^ Z* ^ ZH — > i^£~bb, where 
£ is e or fj,, with a dataset of 1 fb^^, almost three times 
that of the previously reported analysis i^]. CDF and 
D0 have previously presented Hi ggs boson searches in 
other decay modes 0, i, S, M, El, M [ll ■ 



CDF II [lj| is a general purpose detector. Its coordi- 
nate system and quantities used throughout this paper 
are defined in Ref. [15]. At its center is a cylindrical sili- 
con detector which tracks charged particles from a radius 
of 1.35 to 29 cm for \r]\ < 2. Around this is a cylindrical 
wire drift chamber which tracks charged particles from 43 
to 132 cm for \r]\ < 1.3. A superconducting solenoid sur- 
rounds the tracking volume providing a 1.4 T magnetic 
field for momenta measurements. Segmented electro- 
magnetic and hadronic sampling calorimeters surround- 
ing the solenoid measure energies of interacting particles 
with |?7| < 3.6. A system of drift chambers and scin- 
tillation counters outside the calorimeters detect niuon 
candidates for |r/| < 1.5. 

At the Tevatron the cross section for qq ^ Z* ^ ZH 
production for a Higgs boson with mass mn — 115 
GeV/c^ is 1.04 pb ^6j], and the branching ratio B{H -^ 
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hh) is 73% fl?']. To identify candidate ZH events, we 
first search for Z candidates decaying to electron or muon 
pairs. The full selection criteria are described in Ref. [l8{ : 
the most salient features are described here. Events are 
collected using a trigger which identifies a primary elec- 
tron (muon) with Et > 18 GeV (pr > 18 GeV/c) within 
the central region \r]\ < 1.0. The requirements for the 
second electron are relaxed to Et > 10 GeV in the cen- 
tral region, and maintained at 18 GeV for electrons with 
1.0 < \r]\ < 2.4. The second muon must have px > 10 
GeV/c. Energy deposits from leptons must be isolated 
from other energy deposits within AR < 0.4. 

From the measured lepton energies and momenta, we 
reconstruct the invariant mass of the Z candidate and 
require it to be between 76 and 106 GeV/c^. This re- 
quirement is 92% efhcient for real ZH candidates, but 
helps remove non-Z backgrounds. We require oppositely 
charged leptons for muon pairs and for electron pairs 
when the second electron has \ri\ <1.0 due to improved 
tracking efficiency in the central region. 

Higgs boson candidates are then selected by requir- 
ing a jet with Et > 25 GeV and an additional jet with 
Et > 15 GeV, both with |?7| < 2.0. Jets are corrected 
for calorimeter response, multiple interactions, and en- 
ergy loss in the uninstrumented detector regions [19]. To 
enhance signal significance, we implement an algorithm 
to identify the decay of a long-lived hadron containing 
a b quark by reconstructing a significantly displaced sec- 
ondary vertex [l^, [2lj. The efficiency is 40 to 50% for 
6-quark jets and 1 to 3% for u, d, s, or g (light parton, 
l.p.) jets. We consider events in which one or both jets 
are "tagged" by this algorithm. 

Backgrounds originating from W , Z , and tt production 
are determined using leading order Monte Carlo (MC) 
calculations, normalized to next to leading order, fol- 
lowed by a detailed simulation of the CDF II detector. 
We model the Z + bb, Z + cc and Z + l.p. processes by 
first producing the exact leading order multiparton final 
states with the ALPGEN 22] MC program, and then using 
the HERWIG [23| MC program to model the hadronization 
and parton showering. In addition we use an inclusive 
PYTHIA [2J] Z MC sample to compare with the observed 
data and evaluate systematic uncertainties. We model 
ZZ, ZW, and ti background contributions using the 
PYTHIA MC program. A "fake lepton" background arises 
from jets being misidentified as leptons, and we estimate 
this contribution from observed data ^25j . The contribu- 
tion from false tags of l.p. jets is evaluated from data 
by applying a parametrization of the false tag ging rate 
to jets passing Et and rj requirements 20|, |2l|. The ac- 
ceptance, with statistical uncertainties, of ZH -^ £^£^bb 
events is (10.8 ± 0.1)% for mn = 120 GeV/c^ and is 
evaluated using the pythia MC program followed by a 
detailed CDF II simulation for Higgs boson masses from 
110 to 150 GeV/c2. 

TableUshows the expected background and signal con- 



TABLE I; Expected and observed numbers of events in 1 fb~^ 
for electron and muon decay modes combined, compared to 
the expected ZH signal for mH ~ 120 GeV/c^. 



Sample Single-tagged Double-tagged 



Z + hb 


35.1 ±14.6 


6.3 ±2.5 


Z + cc 


21.8 ±8.5 


1.0 ±0.4 


Z+l.p. 


32.3 ±5.5 


1.0 ±0.2 


tt 


5.2 ±1.0 


2.8 ±0.6 


zz 


4.0 ±0.8 


1.3 ±0.3 


zw 


1.2 ±0.2 


0.04 ±0.01 


Non-Z 


1.9 ±1.4 


0.2 ±0.2 



Expected 


101.5 ±32.0 


12.7 ±4.1 


Observed 


100 


11 


ZH 


0.44 


0.23 



tributions with systematic uncertainties in 1 fb^^ of data 
compared to the number observed after dividing events 
into those with only one 6-tag (single-tagged), and those 
with exactly two &-tags (double-tagged). Electron de- 
cay modes account for 60% of the total expected and 
observed events. 

The Z+jets system is only expected to have ^t due 
to jet mismeasurement, either from the limited calorime- 
ter resolution, uninstrumented regions, or from the semi- 
leptonic decay of the jets. Therefore, we correct the jet 
energies using a multilayer perceptron Artificial Neural 
Network (ANN) [26] which uses the missing transverse 
energy ^t vector projected onto those of the jets in or- 
der to determine individual scale factors for each jet. The 
result is a dijet mass resolution improvement from 18% 
to 11% 18] which improves ZH signal discrimination 
from the backgrounds. 

To achieve a greater separation of signal and back- 
ground we employ an additional ANN implemented with 
JETNET 27| to distinguish the kinematics of the signal 



from those of the backgrounds. Our ANN configuration 
is 8 input variables, 17 hidden nodes, and 2 output nodes, 
such that the output distribution is two-dimensional (2- 
D), with one axis separating ZH and Z ± 66, and the 
other axis separating ti and ZH . The variables chosen, 
in order of importance for minimizing the classification 
error, are the scalar sum of the transverse energies of the 
jets and leptons composing the Higgs and Z candidates 
(Ht), ^t, dijet mass, AR between first jet and Z candi- 
date, AR between subleading jet and Z candidate, AR 
between leading and subleading jets, sphericity, which is 
a measure of how isotropic the leptons and jets are, and i] 
of the subleading jet. The most important distributions 
are shown in Fig. [TJ 

The uncertainty on the amount oi Z + bb and Z + cc 
background is taken to be 40% f25|| and for the tt, WZ, 
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FIG. 1: Expected and observed distributions for the three 
most important inputs to the signal discriminating ANN 
shown after ANN jet corrections have been applied, for events 
with 1 6-tag. 



and ZZ [28|, |29|, it is 20%, including the uncertainties 
on the cross section and on the selection efficiencies of 
these processes and the top quark mass uncertainty. The 
uncertainty on the non-Z background is 50%. The un- 
certainty on the shape of the background is evaluated by 
comparing Z + bb events between PYTHIA and ALPGEN. 
The signal shape uncertainty is evaluated by varying the 
amount of initial and final state QCD radiation [30|, and 
by changing the parton distribution functions using the 
40 eigenvectors from CTEQ6 [3l|. We evalute both rate 
and shape uncertainties for the signal and backgrounds 
by varying the jet energy scale within its uncertainties 
[19J . In addition both signal and background estimates 
are affected by the trigger efficiency uncertainty (1%), 
and the luminosity measurement uncertainty (6%) [32] ■ 
The 6-tagging efficiency has an uncertainty of 8% for b- 
quark jets, 16% for c-quark jets and 13% for l.p. jets. 
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FIG. 2: Expected and observed distributions for the ANN dis- 
criminant projected onto the axis which discriminates Z+jets 
from ZH, after enriching signal by selecting the most signal- 
like 25% of events as determined by the ii vs. ZH ANN out- 
put. The top plot is the ANN discriminant before 6-tagging 
requirements, demonstrating two models for the Z-j-jets back- 
ground normalized to data for shape comparison. The mid- 
dle and lower plots are for single-tagged and double-tagged 
events, and are shown with a ZH signal at the level of the 
observed 95% CL. upper limit for itih = 120 GeV/c^. 



In double-tagged events, the uncertainties are 16%, 32% 
and 24%, respectively. These values are updates obtained 
from the procedure found in Ref . [20[ . 

The projections of the 2-D ANN signal discriminant 
are shown in Fig. [2] We analyze the binned 2-D ANN dis- 
criminant distribution to test for a ZH signal in the pres- 
ence of SM backgrounds using a Bayesian technique [33| 
and marginalize over variations in the systematic uncer- 
tainties. The expected and observed upper limits at the 
95% CL. are shown in Tabic [HI Expected limits are 
obtained by generating pseudo-experiments from the ex- 
pected SM ANN shapes to calculate the median ZH con- 
tribution which could be excluded at the 95% level with 
no ZH signal present. Since backgrounds are smaller for 
double-tagged events, we analyze them separately from 



TABLE IL Upper limits at the 95% CL. on the cross-section 
of a{ZH) X B{H -^ bb). Also shown is the expected and 
observed ratio of the limit compared to the SM cross section. 
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single-tagged events, resulting in a 20% improvement in 
the expected limits. 

The dominant systematic uncertainty is the 6-quark 
identification efficiency which accounts for 12% of the 
total 14% increase in the expected limit due to systematic 
uncertainties. 

In summary, we have extended the limits for a Higgs 
boson decaying to hh produced in association with a Z 
boson. This is the first Tevatron Run II search for a Higgs 
boson to use a multivariate approach to separate signal 
and background kinematics, and results in a significant 
improvement in Higgs sensitivity over previous analyses 
in this decay mode [y|. The observed event kinematics 
and ANN signal discriminants show no significant excess 
above SM predictions. The improvement in limits using 
our approach of two ANNs is a factor of 1.8 compared 
to a fit of the uncorrected dijet mass distribution alone. 
This result finds the best limit on Standard Model Higgs 
production in the most favored Higgs boson mass range 
of less than 125 GeV/c^ that was achieved with 1 fb^^ of 
data as in Ref. [l8|. Full sensitivity for this mass range 
will be achieved by combining this analysis with other 
CDF and D0 Higgs search channels and the combined 
Tevatron dataset. 
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